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METHOD AND APPARATUS FOR QUANTIZING A COLOR IMAGE 

THROUGH A SINGLE DITHER MATRIX 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application incorporates by reference to and priority claimed from U.S. Provisional 
Patent Application No. 60/179,586 entitled "Quantizing Color Image Through a Single Dither 
Matrix," filed February 1, 2000. 

BACKGROUND OF THE INVENTION 

Field of Invention 

This invention relates generally to imaging processing and more particularly to color 
image quantization technology. 

Description of the Prior Art 

The present invention concerns the technology of image output to a video display or to a 
printer. In most cases, a digitally stored color image is stored in terms of nominal values that a 
device outputting the image should attempt to approximate as closely to the original source as 
possible. The image may be stored as pixels wherein the pixels are typically presented in the 
form of a three-component vectors. Quantizing or dithering an image is a process of reducing 
bit-resolution of a source image. Image quantization or dithering can be done through either 
error-diffusion or dither matrix, which is also known as the threshold array method. 

Conventionally, quantizing or dithering a red-green-blue (RGB) or cyan-magenta-yellow 
(CMY) color image is done by quantizing each colorant separately through multiple independent 
dither matrices. The method is simple to implement, but at a cost of local color variation and the 
complexity of providing multiple independent dither matrixes. 

Referring to Fig. 1, which illustrates a prior art schematic view of a source color space 10 
of a color image. The color space 10 is presented as a three-dimensional vector space defined by 
three axes A 12, B 14 and C 16 which correspond to colorants R, G, B or C, M, Y respectively 
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depending on the output device. A pixel value is a three-component vector having a point in the 
color space. A color image is a two-dimensional array of pixel values. A source image contains a 
source color image. An output image for display or printout contains an output color image. 
Each pixel in the source image is referred to by Pixel[ij], wherein i and j are indices into the 
5 pixel array. The pixel values of colorants in the source image are referred to as A s [i,j] 1 8, B s [i,j] 
20, C s [i j] 22 respectively. The values of each colorant are integers encoded in S-bits, which 

means 0 <= A s [i,j], B s [i j], C,[i j] < 2 s . 

Fig. 2A, 2B and 2C are diagrams showing a prior art method for dithering a source 
image. The color image is to be dithered into a 3 -colorant pixel array with colorant values 
1 o encoded in T-bits, wherein T <= S. Prior to dithering the source image, the S-bits source 

colorant values A s [i,j], B s [i j] and C s [i,j] are converted to higher bit-resolution encoded in U-bits, 
o where U>=S, through three linearization look-up tables L a 24, L b 25 and L c 26 for one or more of 
Cj the following purposes: linearization, gray balance, and preserving numerical precision. The 
/:? converted U-bits colorant values are in a quantization space (0, 2 U ). The converted three colorant 

H " " 

Gl 5 values in the quantization space are denoted by A q [i j] 27, B q [i ,j] 28 and C q [i j] 29. As shown in 
5 Figs. 2 A, 2B and 2C, three quantization table A t [t,] 30, B t [t 2 ] 31 and C t [t 3 ] 32 for each colorant A, 
B, C respectively are defined, where 0 <= t l5 1 2 ,t 3 < 2 T . The values of A^], B t [t 2 ], and C t [t 3 ] are 

H 1 U-bit numbers in the quantization space. This means 0 <= A^tJ, B t [t 2 ], C t [t 3 ] < 2 , and A t [tJ, 

bl 4. r 

p B t [t 2 ], C t [t 3 ] are selected from the 2 U possible values and satisfying A^O] = 0, B t [0] = 0, C,[0] = 0; 
H20 AJ.t^A.ft+l], B t [t 2 ]<B t [t 2 +l], C t [t 3 ]<C t [t 3 +l]. It is desirable, but not required, to have a equally 

divided quantization space, that is for all the possible values of a and b, (A t [a+1]- 
A t [a])KA t [b+l]-A t [b]) ? (B t [a+l]-B t [a])=(B t [b+l]-B t [b]) and (C t [a+l]-C t [a])KC t [b+l]-C t [b]). The 
indices t x 36, t 2 37 and t 3 38 which are indices of the three quantization tables are also the output 
pixel value in the output color space. 

25 Fig. 3 is a diagram showing a prior art schematic view of a 3 -dimensional cubic 

quantization interval 40 containing the converted pixel value (A q [i j], B q [i j], C q [ij]) 41 in the 
quantization space 42. The entire quantization space are divided into smaller cubic subintervals 
by the three quantization tables A t [ ], B t [ ] and C t [ ]. As shown in Fig. 4, which is a diagram 
illustrating a prior art schematic view of a cubic quantization interval, having 8 corner points 

30 P 0 oo=A t [t 1 ]B t [t 2 ]C t [t 3 ] 44, P 100 =A t [t 1 +l]B t [t 2 ]C t [t 3 ] 46, P 010 =A t [t I ]B t [t 2 +l]C t [t 3 ] 48, 
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Pno=A t [t,+l]B t [t 2 +l]C t [t 3 ] 50, P fl01 =A t [t 1 ]B t [t 2 ]C l [t,+l] 52,P 101 =A t [t 1 +l]B t [yC t [t 3 +l] 54, 
Po,rA t [t,]B t [t 2 +l]C t [t 3 +l] 56, P„,=A t [t 1 +lMt 2 +l]C t [t 3 +l] 58. 

A dither matrix, or threshold array is an two dimensional array of uniformly distributed 
5 random or psudo-random numbers H[i j], where 0<=H<=H max . If the quantization space is 
equally divided, equals to (A^l]-^,]), (B^l]^]) and (C t [t 3 +1]-C t [t 3 ]). 

For each pixel in the color image, the goal of the dithering process is to select one out of 
the eight corner points P 000 , P ]00 ,. . ., P,„ of the cubic quantization interval 40 containing the 
10 converted pixel value 41 in the quantization color space 42, by comparing a value from the 
threshold array against the relative offset of the pixel value in the quantization interval. The 
n output pixel value of the dithering process is three quantization table indices t„ t 2 and t 3 
i;1 corresponding to the selected comer point. 

Q5 The above-described conventional method of quantizing or dithering an image requires 

rjj three independent dither matrices for each of the three colorants to process the source image. 
!L The result of this methodology tends to generate color with higher variation. It is, therefore, 
H desirable to have a method that can efficiently dither an image with a single dither matrix yet still 

i* * 1* 
»> **■ 

£1 be able to retain the quality of such image. 

[720 

SUMMARY OF THE INVENTION 

An object of the invention is to dither a color image with smaller color variation with a 

single dither matrix. 

Another object of the invention is to dither a color image with better control of dot 
25 placement of multiple colorants. 

In a specific embodiment of the presentation, a device for dithering a color image is 
described. In one embodiment, the device includes a truncating logic for truncating a value of a 
pixel having a point in a color space of said source image, for locating a 3 -dimensional cubic 
subinterval containing the point; a random number generator for generating a random number as 
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a threshold value for dithering each pixel value of said source image; a first comparator for 
locating a tetrahedron having four vertices in said three-dimensional cubic subinterval that 
contains the pixel point; a second logic for assigning a probability density corresponding to each 
of the vertex; a second comparator for selecting a vertex provided by the threshold value through 
an inverse probability distribution function derived from the probability densities of the vertices; 
and a third logic for locating an output point in an output color space corresponding to the 
selected vertex. 

In another embodiment, a device for dithering a source color image via weighting 
coefficient associated with a plurality of sample points in a source color space of the source color 
image is described. The device includes a first logic for generating the weighting coefficients of 
the sample points, wherein the weighting coefficients are probability densities of the sample 
points; a second logic for generating a random number as a threshold value for dithering a pixel 
value of the source image; a third logic selecting a point from the plurality of said sample points 
provided by the threshold value through an inverse probability distribution function derived from 
the weighting coefficients of said points; and a fourth logic for locating an output point in an 
output color space corresponding to the selected point. 

One advantage of the present invention is that it provides a mechanism that is capable of 
dithering a color image with smaller color variation or error. 

Another advantage of the present invention is that it provides a mechanism that is capable 
of dithering a color image with better control of dot placement of multiple colorants. 

These and other objects and advantages of the present invention will no doubt become 
apparent after reading the following detailed description of the preferred embodiments which are 
illustrated in the several figures of the drawing. 

TN THE DRAWINGS 

Fig. 1 is a diagram illustrating a prior art schematic of a color space of a source color 

image; 

Figs. 2A-2C are diagrams depicting a prior art method of dithering a source image; 
Fig. 3 is a diagram showing a prior art schematic view of a pixel point in a cubic 
subinterval of a quantization space; 
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Fig. 4 is a diagram illustrating a prior art schematic view of a cube subinterval of a 

quantization space; 

Fig. 5 is a block diagram illustrating an exemplary embodiment for the image processing 
method and apparatus of the present invention; 

Fig. 6 is a block diagram illustrating the software operation and flow in the embodiment 

depicted in Fig. 5; 

Fig. 7 A is a diagram illustrating a method of one-dimensional interpolation using 
difference values; 

Fig. 7B is a diagram illustrating a method of a two-dimensional interpolation using 

difference values; 

Fig. 7C is a diagram showing a three-dimensional interpolation unit; 

Figs. 8A and 8B are diagrams showing two forms used in triangular interpolation in a 2- 

dimensional color space; 

Fig. 9A-9C are diagrams each illustrating a different way of dissecting a 3-dimensaional 
cubic subinterval into tetrahedrons and showing weighting coefficients associated with the 
vertices; 

Fig. 9D is a diagram illustrating a lattice structure of tetrahedrons; 
Fig. 10 is a diagram illustrating a method of deriving weighting coefficients of a 
tetrahedron of the preferred embodiment of the present invention; 

Fig. 1 1 is a block diagram illustrating an exemplary dithering process of the present 

invention; 

Fig. 12 is a diagram showing schematic view of a pixel value in a cubic subinterval in a 

color space of the present invention; 

Figs. 13 A and 13B are diagrams illustrating the probability density function and 
associated probability distribution function used in the present invention; 

Fig. 14 is a block diagram showing an exemplary embodiment of the structure for the 

present invention; and 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The present invention discloses a system and method for generating a color image for 

outputting to an output device such as a video display, printer, facsimile machine, etc. The 
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system according to the invention converts a source image into an image that has less pixel depth 
and yet preserves the visual image information of the original source image. The invention can 
be embodied either in a computer program that implements the procedure described later in this 
invention, or in dedicated digital electronic circuitry designed to perform the same function. 
5 The present invention utilizes a statistical algorithm in converting pixel values in the 

source image into output values through quantizing RGB or CMY colorant by selecting an 
corner point of a region from a plurality of predefined regions in color space that is closest to the 
pixel to be the output of the pixel of the image. 

Fig. 5 is a block diagram illustrating an exemplary embodiment for the image processing 
10 method and apparatus of the present invention. A computer system 60 sends formatted image 
data stream to a multi-color display devices such as a CRT display 62 or color printer 64 via 
O appropriate communications channels such as a video cable 66 and a printer cable 68 
CI respectively. The printer cable 68 can be a parallel communications cable or a serial 
^ communications cable, or a USB communications cable. The computer system 60 comprises 
Ol5 various devices for generating the formatted image stream such as a central processing unit 
hi (CPU) 70 and memory 72 coupled together through a bus 74 in the computer system 60. The 
L. computer system 60 further includes a display adapter 76 coupled to the bus 74 on one end and to 

II I . I 1 1. 

N= the CRT display 62 via the video cable 66. The CPU 70 processes the source image and 

I-" -"A M 

********** 4 _^ . 

£1 generates the output data to be displayed on the CRT display 62 or to be printed on the printer 64 
fJ20 in accordance with the algorithm to be discussed hereunder. The computer system 60 may 

further include an input/output interface (I/O interface) 78 for handling communications between 
the computer system 60 and the printer 64, or a facsimile machine (not shown). 

In the embodiment described above, the printer 64 comprises a processor 80, a memory 
82, and a control logic 84 to control the printer operation as well as the communications to and 
25 from the I/O interface 78. The processor 80 is coupled to the memory 82 and the control logic 84 
and controls the operation of the control logic 84. The processor 80 may further process the 
source image it receives from the computer system 60 and operative to produce the output data to 
be printed on the printer 64 according to the algorithm to be discussed below. 

Fig. 6 is a block diagram illustrating the software operation and flow in the embodiment 
30 depicted in Fig. 5. The computer system 60 comprises an application software 86 that performs 
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the operation in converting a source image to an output data to be displayed on a video display or 
printed on a printer. The application software sends the converted output data to a buffer where 
display driver 88 in the display adapter 76 can retrieve for outputting to the CRT display 62. 
Alternatively, the application software 86 may send the converted output data to a buffer where a 
printer driver 90 embedded in the I/O interface 78 may retrieve for printing to a printer 64. The 
printer 64 includes an embedded software 92 that controls the operation of the control logic 84 in 
the printer for printing the output data on the printer 64. In the preferred embodiment, either the 
CPU 70 in the computer system 60 or the processor 80 in the printer can employ the 
methodology described in the present invention to convert pixilated source image data to lower 
pixel -depth suitable for the destination output device. As depicted in Fig. 6, in a typical 
implementation the instructions performing the invention's method are usually embedded in a 
program such as display driver 88 or printer driver 90, or as part of an embedded printer software 
92. 

For better performance concerns, the display adapter 76 or the printer control logic 84 can 
implement present invention's methodology in a digital logic circuitry which can perform the 
operation in much higher speed. In such implementation the CPU 70 of the computer system 60 
sends the unconverted source image to the image displaying device for conversion. 

As a background to the present invention, the method of linear interpolation will be 
described. A set of probability density functions described later in the invention is derived from 
the general linear interpolation function, which is: 

(V) = I a, * (V,) 

* 

i 

I a ; = 1 and 0 < a, < 1 

1 

where 
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"Vi 's" are the vertices that defines a subinterval area over which 
the interpolation is performed. Each V; has a corresponding output value 

(VD. 

"V" is the input point for which the interpolated output value (V) is 
derived. V is located in the subinterval area defined by Vj 's. 

"a," is the weighting coefficient corresponding to the vertex V { . a t 
is directly derived from the relative position of V with respect to the 
positions of the vertices V r a L meets the requirements that the closer V is 
located to a particular vertex V x , the value of a x becomes larger. If V 
happens to be one of the vertex V x , the value of a,, is 1 and the values of all 

the other a, (x * i) are 0. 
Referring now to Fig. 7 A, which is a diagram showing a one-dimensional interpolation 
method. In Fig. 7 A, V 0 94 and Vj 96 denote two input values which defines a subinterval of a 
one-dimensional space. Each of these input values V 0 94 and V\ 96 maps to a corresponding 
output value (V 0 ) 98 and (Vj) 100 as shown. Assuming the distance between V 0 and V, is 
normalized to 1, the interpolation for any input value that falls between V 0 and V\ can be 
expressed with a weighting coefficient "a" 102, 0 <= a <= 1, wherein "a" is the normalized 
distance between V and V 0 . The normalized distance between V and V 5 is (1-a) 104. Based on 
the location of a input value V that falls between V 0 and V 1? an interpolated output value (V) 
corresponding to the input value V located in the subinterval between V 0 and V, is given as 
follows: 

(V) = (l-a)*(V 0 ) + a*(V 1 ) 
where a is the normalized distance from V 0 to V: 

a - (V - V 0 )/( V 1 - V 0 ) 
The basic interpolation form can be extended into multi-dimensional space. Fig. 7B is a 
diagram depicting a rectangular area subinterval in the 2-dimensional space. The input values at 
the four corners are labeled V 00 106, V 10 108, V 01 1 10 and V n 1 12 with the corresponding output 
values are (VJ 114, (VJ 116, (VJ 118 and (VJ 120. The interpolated output value (V) 122 
of any input value V 124 in the 2-dimensional subinterval defined by V 00 , V 01 , V 10 , V n is given 
by using two weighting coefficients "a" 126 and "b" 128 as follows: 
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(V) = f (x 5 y) - (1-a) * (1-b) * (V 00 ) + a * (1-b) * (V 10 ) + 

(l-a)*b*(V 01 )+a*b*(V u ) 
where a 126 and b 128 are the normalized distance from V 00 to V: 

a = (V x -V 00iX )/(V Rx -V 00)X ) 
b-(V y -V 005y )/(V oliy -V 00)y ) 

where V x is the position of V in reference to the X-axis and V y is the 

position of V in reference to the Y-axis. 
Referring now to Fig. 7C, which is a diagram showing a three-dimensional subinterval 
over which interpolation is performed. Similar to the above description in 1 -dimensional or 2- 
dimensional interpolation, in a 3 -dimensional space, over the cubic subinterval defined by 
vertices V 000 126, V 100 128, V 010 130, V 110 132, V 00l 134, V 101 136, V 011 138, V m 140 the 
interpolated output value (V) (not shown) of an input value V 142 in the subinterval has the 
following form: 

(V) = (1 - a) * (1- b) * (1 - c) * (V 000 ) + a Ml - b) * (1- c) * (V 100 ) + 

(1- a) * b * (1 - c) * (V 010 ) + a * b * (1- c) * (V ]10 ) + 
(1- a) * (1 - b) * c * (V 001 ) + a * (1-b) * c * (V 101 ) + 
(1- a) * b * c * (V ou ) + a * b * c * (V m ) 
where a, b, c are normalized distance from V 000 to V: 
a = (V x - V 000i x ) / (V 100 , x ~~ ^ooo, J 
b = (V y - V 000i y ) / (V 010f y - V 00 o 5 y ) 
c = (V 2 - V 000i z ) / (V 001) z - V 000i J 

where V x is the position of V in reference to the X-axis, V y is the position 
of V in reference to the Y-axis, and V z is the position of V in reference to 
the Z axis. 

To reduce the numerical operation involved in interpolation process, it is common 
practice to use less sample points or vertices in the process, yet still obtains acceptable results. In 
the 2-dimensional space, a minimum of three sample points that define a triangular shape 
containing the interpolation point V are required to produce the interpolated output. In the 3- 
dimensional space, a minimum of four sample points that define a tetrahedron volume containing 
the interpolation point V are required to produce the 3 -dimensional interpolated output. 
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In a 2-dimensional space as illustrated above, there are two forms of triangular 
interpolation for computing the interpolated output. One is shown in Fig. 8A, the interpolated 
output value (V) is given by arranging weighting coefficient for the upper triangle 144 and lower 

triangle 146 as following: 

(V) - (1-a) * (V 00 )+ (a-b) * (V 10 ) + b * (V„) if a > b; or 

(V) = (1-b) * (V 00 ) + (b-a) * (V 01 ) + a * (V„) if a < b 

Referring to Fig. 8B, in the alternative interpolation form the 2-dimensional rectangle 
subinterval is dissected into upper triangle 148 and lower triangle 150 by drawing a line through 
vertices V 10 and V 01 . The interpolated output is given as following: 

(V) = (1 - a - b) * (Voo) + a * (V 10 ) + b * (V 01 ) if (a + b) < 1; or 

(V)= (l-b)MV ]O ) + (l-a)*(V 0 0 + (a + b-l) ;i; (V 11 ) if(a + b)>l 
The tetrahedral interpolation forms use 4 corner vertices out of the 8 vertices of a cubic 
subinterval in the 3-dimensional space. There are many 3-dimensional tetrahedral interpolation 
forms. Figs. 9 A, 9B, and 9C are diagrams each illustrating one of 3 different ways of dissecting 
a 3-dimensional cubic subinterval into tetrahedrons. A list of the weighting coefficient of the 
vertex subsets of each tetrahedron is also given for each interpolation form. Fig 9D is a diagram 
illustrating a lattice structure of tetrahedrons which gives a more uniform tetrahedrons but 
requires more complicated interpolation coefficient calculations. 

The manner in which the four weighting coefficient a, of a tetrahedron is derived as 

follows. 

Fig. 10 is a diagram illustrating a method of deriving weighting coefficients of a 
tetrahedron of the preferred embodiment of the present invention. The coefficient is the 
probability for the output data to fall within the tetrahedron within a dissected unit cube. 
Referring now to the example illustrated in Fig. 10, first select one of the four vertex as the new 
origin point W 152 identified as V 100 154, and three axis X 156, Y 158, Z 160 pointing to the 
other three vertices identified as V on 162, V 010 164, and V 001 166 respectively. Next, write down 
the vector representation of W, X, Y, Z in terms of the base vectors A 168, B 170, C 172 and the 
origin O 174. Follow by representing O, A, B, C in terms of W, X, Y, Z. The last step is to 
convert the representation of an input interpolation point from the O, A, B, C space into the new 
W, X, Y and Z space. The coefficients of X, Y, Z in the new representation are the weighting 
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coefficient of vertices X, Y, Z respectively. The weighting coefficient of W is given by 
subtracting the weighting coefficients of X, Y, Z from 1. 

In the examples of Fig. 10, W, X, Y, Z are written as: 

W = 0 + A 
X = C-A 
Y = B - A 
Z=B+C-A 
Then solving for 0, A, B, C: 
0 = W-A 
A=Z-X-Y 
B = Z-X 
C = Z-Y 

Substitute O, A, B, C into the vector representation of an interpolation point: 
O+a* A+b*B+c*C 

= (W- A) + a*A + b*B + c*C 

- W + (a-l)* A+b*B+c*C 

= W + (a-l)*(Z-X-Y) + b*(Z-X) + c*(Z-Y) 

= W + (a + b + c-l)*Z + (l-a-b)*X + (l-a-c)*Y 

The resulting weighting coefficients of vertices X, Y, Z are (1 - a - b), (1 - a - c) and (a 
+ b + c - 1 ) respectively. The weighting coefficient of W is 1 -(1 - a-b) - (1 - a- c) - (a + b + 

c - 1) = a. 

It is the basic feature of the present invention that the 2-dimensional triangular or 3- 
dimensional tetrahedral interpolation coefficients described above are used as probability density 

values for dithering over the tetrahedron within which the input point is located. 

Fig. 1 1 is a block diagram illustrating an exemplary dithering process of the present 
invention. Input color pixels are encoded as 3-colorant values. Each colorant is encoded in an S- 
bit binary number, wherein the S-bit is the original precision of the colorant value. The input 
color image is to be converted into an output image of lower bit depth T for each colorant. The 
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colorants are named A, B and C where in a particular application A, B, C can refer to R, G, B 
display color or C, M, Y printer color, or any device independent color such as L, a, b color. An 
image source provides a 3-colorant source image 176 to be converted into a dithered output 
image 178. The image source can be a scanner, digital camera, or image data stored in computer 
memory. Three 1 -dimensional look-up tables 180 each converts its respective S-bit input 
colorant value to a higher bit depth U, where U > S. This process avoids the loss of precision 
before the quantization occurred in conventional color calibration processes. In conventional 
color calibration processes, an 8-bit to 8-bit look-up table can map two input color values onto 
the same output value and, thus, causes lost of precision as a consequence. These look-up tables 
are optional, but preferred in a typical implementation of the present invention for one or more of 

the following practical reasons: 

1 . normalizing the quantization intervals; 

2. linearizing or calibrating output device; 

3. maintaining gray balance from the input to the output image; or 

4. preserving numerical precision during subsequent computation. 

The first reason identified above eliminates the need to build 3 quantization tables of non- 
uniform quantization intervals. These tables of non-uniform quantization interval are used to 
locate the 3 -dimensional cubic subinterval that contains the input pixel value. With normalized 
quantization intervals the cubic subinterval can simply be located by truncating off the lower bits 
of the U-bit colorant values and use the higher bits as the indices to the cubic subinterval. 

Still referring to Fig. 1 1, the truncation functions 182 separates the upper T-bits of the U- 
bits color values and the remaining lower order (U-T) bits. The upper T-bits of the 3 colorants 
identify the cubic subinterval that contains the input point to be described in conjunction with 
Fig. 12. As illustrated in Fig. 7C and Fig. 12, the value of each pixel corresponds to a point 41 in 
the 3-dimensional color space 42. The vertices of the cubic subinterval containing the pixel point 
are named V 000 , V 100 , V 010 , V n0 , V 00lJ V 101 , V 010 , V 1U , 

The three lower-order (U-T) bit values are named a, b, c for each colorant A, B, and C 
respectively. As illustrated in Fig. 7C in conjunction with Fig. 12, a 193, b 194 and c 195 are the 
relative position of the pixel point 41 in the cubic subinterval 40 in the color space 42 with 
respect to the origin of the lower corner V 000 196 of the cubic subinterval. The lower-order (U-T) 
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bits of a U-bit colorant value represent the relative position of the pixel value within the cubic 
subinterval identified by the upper T-bit value. Since the value of a, b 5 c are (U-T) bits, the unity 
value 1 representing the size of the cubic subinterval is scaled up by 2 (U " T} times and denoted by 
IN for the subsequent integer computation. 

Referring to Fig. 1 1, A pseudo-random number source 184 provides a uniformly 
distributed threshold value H, where 0 < H < 2 (U ~ T) . The pseudo-random number source is usually 
implemented as a dither matrix, or threshold array, which is tiled and overlaying on top of the 
source image. The values of the dither matrix are periodically fetched from the matrix when 
each pixel of the source image is fetched and processed. In the case of bi-level output image 
values (T = 1), the types of dither matrix are classified in two categories for practical reasons. 
The clustered-dot dither matrices tend to have high element values clustered together. The 
dispersed-dot dither matrices tend to have high-valued matrix elements mixed with low-valued 
ones. Dispersed-dot dither matrices in general produce better output quality over the clustered- 
dot dither matrices, which perform better when the output device has a very nonlinear response 
with respect to the input pixel values. 

A set of comparing functions 186 identifies the tetrahedron in the identified cubic 
subinterval further contains the input point. The actual comparing functions vary in accordance 
to the chosen tetrahedral interpolation form. Given the tetrahedral interpolation form illustrated 
in Fig. 9C, as an example, the following numeric comparisons are performed: 



Table 1. Comparing functions to determine a tetrahedron for Fig. 9C. 



Comparisons performed 


Tetrahedron/vertices identified 


(a + b + c) < IN 


Ti / V 100 , V 010 , V 001? V o00 


(a + b + c) > IN, (a + b) < IN, (a + c) < IN 


T2 / V 100? V 0]0? V 00] j V 011 


(a + c) > IN, (a + b) < IN 


T3 / V 100 , V J01 , V 001 , V 01l 


(a + b) > IN, (a + c) < IN 


T4 / V 100 , V 010 , V 110? V 011 


(a + b + c) < (2* IN), (a + b) > IN, (a + c) > IN 


T5 / V 100 , V 101 , V n0 , V 011 


(a + b + c)>(2*lN) 


Te / V 1U? V 101 , V il0 , V 011 
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The tetrahedron identified by the comparing functions further determines the set of 
weighting coefficients used in the dithering. Referring to Fig. 1 1, a weighting coefficient table 
188 is constructed by rearranging the weighting coefficients of the chosen interpolation form. 
Using the same interpolation form as illustrated in Fig. 9C, the weighting coefficient is arranged 
in table 2 as following: 



Table 2. Listing of the vertices of each tetrahedron and the 
corresponding weighting coefficients as shown in Fig. 9C. 



Tetrahedron 


v,/p, 


V 2 /P 2 


V3/P3 


V 4 /P 4 


T, 


v 100 


a 


v 010 


b 


v 001 


c 


v 

v 000 


1-a-b-c 


T 2 


V.oo 


a 


V 0 ,o 


1-a-c 


v 001 


1-a-b 


V 011 


a+b+c-1 




v 100 


1-c 


v 10I 


a+c-1 


V 001 


1-a-b 


V 01I 


b 


T4 


v 100 


1-b 


V 0I0 


1-a-c 


V„ 0 


a+b-1 


Von 


c 


T 


V 100 


2-a-b-c 


v 101 


a+c-1 


V„ 0 


a+b-1 


V 011 


1-a 


T 


v m 


a+b+c-2 


V I01 


1-b 


V„ 0 


1-c 


v 011 


1-a 



For each row of Table 2 ? the actual ordering of the vertices and weighting coefficients 
bears no mathematical difference. However, care must be taken, especially in the clustered-dot 
dither matrix situation, when arranging the weighting coefficients across any two adjacent 
tetrahedrons that share two or three common vertices. Drastic dot pattern change can be avoided 
when smoothly shaded pixel values in the source image transition from one tetrahedron into a 
neighboring one. 

Referring now back to Fig. 1 1 ? the thresholding functions 190 comprises a second set of 
comparators that compare the threshold value H against the weighting coefficients from the 
weighing coefficient Table 2 in association with the selecting and computing the weighting 
coefficient 188. Referring to the probability distribution function depicted in Fig. 13A, the 
weighting coefficients are probability densities of the corresponding vertices V l3 V 2? V 3 and V 4 as 
shown in Fig. 13B. The output random variable V x takes the value of V 1? V 2> V 3 or V 4 and is 
statistically determined by the uniformly distributed random variable H through the inverse 
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function of the probability distribution function. With the weighting coefficient described in 
Table 2, the comparisons with the thresholding functions proceed as follows: 

• First compare H with the first weighting coefficient P t of the chosen tetrahedron (or row) 
from Table 2. If H is smaller than or equal to the weighting coefficient P l5 stop here, the vertex 
V! of the chosen tetrahedron corresponding to the weighting coefficient P, is chosen as the 
output. If H is greater than P ls subtract P, from H, compare the remaining H against the second 
weighting coefficient P 2 in the row, if H smaller than or equal to P 2 , stop here and output vertex 
V 2 of the chosen tetrahedron, otherwise subtract P 2 from H and compare the remaining H with 
the third weighting coefficient P 3 of the chosen row. Again if H is smaller than or equal to P 3 , 
output vertex V 3 of the chosen tetrahedron and stop, otherwise V 4 of the chosen tetrahedron is 
selected. 

When implemented with an electronic circuitry, the comparisons of the thresholding 
functions can be performed in parallel by executing the following compares simultaneously: 

H<P, 

H < P, + P 2 

H<P,+P 2 +P 3 

The resulting vertex within the 3 -dimensional cubic sub interval containing the vertex is 
the dithered output from the dithering procedure. Referring to Fig. 1 1, the last output stage 192 
combines the identified cubic subinterval from the truncation function and the selected vertex 
from the thresholding functions, and gives three T-bits binary values for each of the three 
colorant A, B ? C respectively as the dithering output values to the output destination 178. 

From the probabilistic point of view, various interpolation forms makes no difference in 
terms of giving the same average output colorant values. However the form does have an impact 
to the distribution of colorants with respect to each other. Experiments show that the 
interpolation form given by Fig. 9 A tends to overlay colorants on top of each other, and the 
interpolation form in Fig. 9B and Fig. 9C tends to scatter the colorants. The selection of 
interpolation forms may have a significant effect to the visual quality of the output image. 

* Fig. 14 is a block diagram showing an exemplary embodiment of the structure for the 
present invention. Three colorants of a source image are provided and they are identified as R 
197, G 198, and B 199 respectively in this example. The RGB described here correspond to the 
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three color image source 176 referred to in Fig. 1 1 above. Each of the source colorant values is 
fed to its respective 8-it to 12-bit look-up table 200, 202, and 204. The three look-up tables 
described here correspond to the look-up table 180 described above in Fig. 1 1 . The results 
obtained from each look-up table is fed into a truncation function 206, 208, 210 for each colorant 
5 R,G, B respectively, wherein the upper bits 212, 214, 216 are separated and the lower-order bits 
218, 220, 222 are fed into a comparator 224 to locate a tetrahedron having four vertices in a 
three-dimensional cubic subinterval that contains the pixel point, and to a calculator 226 for 
calculating of weighting coefficients and assigning a probability density corresponding to each of 
the vertices. The comparator 224 and calculator 226 described here correspond to the comparing 
10 functions 186 and computing function 188 described above in Fig. 11 respectively. A random 
number generator 228 having an X and Y position register 230 and a threshold value array 232 is 
r*s provided for generating a random number as a threshold value for dithering each pixel value of 
' yy the source image. The random number generator 228 corresponds to the random number source 
Nl 1 84 in Fig. 1 1 . The threshold value obtained from the random number generator 228 and the 
nl5 output produced by the comparator 224 and the calculator 226 are fed to a comparator unit 234 
H that compares the weighting coefficients to the threshold value for the located tetrahedron and is 
5 operative to select a vertex provided by the threshold value through an inverse probability 

h ■ >ps*i 

ii function derived from the probability densities of the vertices. Finally, the selected vertex is 
rf provided to a logic 236 and combined with the upper bits 212, 214, 216 for locating an output 
D 20 point in an output color space corresponding to the selected vertex and thus determining the 

dithered output value of the source image. The comparator unit 234 and logic 236 described here 
correspond to the thresholding functions 190 and determining logic 192 described above in Fig. 
1 1 respectively. 

While the invention has been particularly shown and described with reference to the 
25 preferred embodiments, it will be understood by those skilled in the art that many other 

modifications and variations may be made thereto without departing from the broader spirit and 
scope of the invention as set forth in the claims. The specification and drawings are accordingly 
to be regarded as illustrative, rather than in a restrictive sense. 

What is claimed is: 
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